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We study the efficiency of a single particle Szilard and Carnot engine. Within a first order 
correction to the quasi-static limit, the work distribution is found to be Gaussian and the correction 
factor to average work and efficiency only depends on the piston speed. The stochastic efficiency is 
studied for both models and the recent findings on efficiency fluctuations are confirmed numerically. 
Special features are revealed in the zero temperature limit. 
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I. INTRODUCTION 

Single particle engines have long been regarded as 
thought experiments, designed to probe the very foun¬ 
dations of thermodynamics and statistical mechanics. A 
celebrated example is the Szilard engine [in , conceived 
to investigate the notion of information in thermodynam¬ 
ics. By now, with the spectacular progress in bio- and 
nanotechnology, these machines are part of an experi¬ 
mentally accessible reality 1M3. Alongside these de¬ 
velopments, there has been significant progress on the 
theoretical side, starting with the discovery of the fluc¬ 
tuation theorem and the Jarzynski and Crooks relations 
[lll - [T^ . and developing into encompassing theories such 
as stochastic thermodynamics [l4l Il5l| . An intriguing re¬ 
cent discovery concerns the properties of the stochastic 
efficiency [13, ll3 - [^ . This new concept derives from the 
observation that the performance, and hence also the effi¬ 
ciency of small machines is subject to strong fluctuations. 
Based upon the fluctuation theorem, it was demonstrated 
that for asymptotically long times, the reversible effi¬ 
ciency is least probable for a machine operating under a 
time-symmetric schedule, while under a time-asymmetric 
driving the probability distributions of efficiency of the 
time-forward and the time-reversed protocol intersect at 
the reversible efficiency. 

In the present paper we investigate a generic set-up 
for a stochastic engine, namely a sin gle p article moving 
between a thermal wall and a piston [26l - l^ . and apply 
the obtained results to the Szilard engine and a simpli¬ 
fied Carnot heat engine. For both engines we investi¬ 
gate the macroscopic efficiency and the efficiency fluc¬ 
tuations. The paper is organised as follows. Section HIl 
introduces the Szilard engine and the simplified Carnot 
engine. While the principle behind each engine is differ¬ 
ent, i.e. information to work conversion and heat to work 
conversion, respectively, the constitutive process, namely 
a particle moving between a stationary thermal wall and 
a perfectly reflecting moving wall, is the same. In sec¬ 
tion Ell we analyse the quasi-static operation including 
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a first order correction. In this limit, the work probabil¬ 
ity distribution for the constitutive process is Gaussian, 
analogous to other isothermal processes [ssl - fs^l . and the 
average work is the free energy difference times a fac¬ 
tor that only depends on the reduced velocity defined as 
the ratio of the velocity of the piston and the thermal 
velocity associated with the thermal wall. The implica¬ 
tions for both engines are discussed. In section IlYl we 
present a detailed numerical analysis of efficiency fluctua¬ 
tions, including finite time results and confirming the pre¬ 
dicted universal properties in the asymptotic time limit. 
A slightly modified heat engine is introduced for which 
analytical results are possible. Conclusions and outlook 
are given in section El 


II. GENERAL SET-UP 

The various heat and work exchange processes are 
schematically represented in figure [TJ We consider, with¬ 
out loss of generality, a one-dimensional setting. The 
constitutive process is shown in figureH^. and corresponds 
to a single particle with mass m moving freely between 
a stationary thermal wall and an infinitely heavy mov¬ 
ing piston. Two types of thermodynamic events occur: 
either the particle collides with the thermalising wall, or 
it collides with the piston. In case of a collision with the 
thermalising wall, the particle gets a new velocity drawn 
from the Rayleigh distribution 

^2 

(j){v)=ve~—, (1) 

with V the reduced, dimensionless velocity o f the par ticle, 
measured in terms of the thermal velocity ^ksT/m with 
kg the Boltzmann constant and T the temperature of the 
wall. The change in kinetic energy then corresponds to 
an exchange of heat with the thermal reservoir. Colli¬ 
sions with the piston are assumed purely elastic. Hence 
a particle with incoming velocity v, hitting the piston 
moving at (dimensionless) velocity u, will recoil with ve¬ 
locity 2u — v. The change in kinetic energy of the particle 
implies an amount of work 2kBTu{v — u) done on the pis¬ 
ton. 
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left compartment, the thermodynamic cost of the mea¬ 
surement is: = —fc^lnx [1^, where x = l/L is 

the compression ratio. The equilibrium distribution is 
now restricted to the left compartment, i.e. given by @ 
but with L replaced by 1. The piston next moves to the 
right with velocity u until the full length L is reproduced, 
and an amount of work w is delivered. For finite piston 
speed, the system will deviate from equilibrium and the 
average of this work is bounded by the decrease AF in 
free energy, {w) < AF with: 

l^F = -kBT\nx (3) 
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Figure 1; Schematic representation of: a) The constitutive 
process: a single particle bounces back and forth between a 
thermal wall on the left and a moving piston on the right, b) 
The Szilard engine: a piston is inserted a distance I from the 
left wall, and subsequently moves away from the particle, c) A 
simplified Carnot heat engine: alternating cycles of expansion 
and compression, while in contact with a hot or cold left wall 
at temperatures Th and F respectively; note the absence of 
an adiabatic expansion and compression phase. 


When the particle is found in the right compartment, 
analogous results hold for the measurement cost Az, the 
decrease in free energy AF and the work w, by replacing 
X with 1 — a;. The efficiency of the above information to 
work conversion is defined as: 


We will investigate the properties of this stochastic quan¬ 
tity in more detail below. Upon repeating the Szilard 
cycle many times, rj will converge to a macroscopic ef¬ 
ficiency fj, which is, as we show below, bounded by the 
reversible efficiency 1: 


ksTiAi)-^’ 


( 5 ) 


see also figure [2^. Referring to the discussion of stochas¬ 
tic efficiency, we also require the time-reversed process 
of the Szilard engine, which is defined as follows. Ini¬ 
tially, the particle is in thermal equilibrium and moving 
between the two thermalising walls. The piston is then 
brought in from one of the sides and pushed inward up 
to position 1. The piston either starts at position 0, with 
probability 1 — a;, or at position L probability x. In both 
cases, work is converted into information. 


B. The simplified Carnot heat engine 


A. The Szilard engine 

For the Szilard engine, shown in figure [lb, the vari¬ 
ous steps per cycle are as follows. At the start, a par¬ 
ticle bounces back and forth between two thermal walls 
at temperature T and separated by a distance L. The 
probability distribution for its position Xq and (dimen¬ 
sionless) velocity vq is given by thermal equilibrium: 

A piston is then inserted a distance I from the left wall, 
trapping the particle either on its left or right hand side. 
A measurement is performed to find out in which com¬ 
partment the particle resides. If the particle is in the 


The various stages of the single particle Carnot heat 
engine are represented in figure [T]:. When discussing the 
heat engine, the reduced velocities are calculated with 
respect to the thermal velocity at temperature Th- The 
cycle starts with the piston moving outward (stage 1) 
at constant velocity u until it reaches its final position 
L. The particle is allowed to relax (stage 2), i.e. the 
equilibrium distribution Eq. is restored. After re¬ 
laxation, the temperature of the heat bath is instanta¬ 
neously switched to the lower value Tc < Th, and the 
particle is again allowed to relax to the new equilibrium 
state (stage 3). Next (stage 4) the piston moves inwards 
with velocity u until it has returned to its original po¬ 
sition I, after which the particle again relaxes (stage 5). 
Finally, the temperature is instantaneously switched back 
to the higher value Th, followed by a hnal relaxation 
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Figure 2: (Color online) Macroscopic efficiency of a) the Szi¬ 
lard engine and b) the heat engine (with = 2Tc) in function 
of u and x. 


(stage 6). The inclusion of the intermediate relaxation 
stages destroys any correlations in position and veloc¬ 
ity that would otherwise appear between the successive 
stages, the relaxation stages are purely dissipative, there¬ 
fore decreasing the efficiency, and have an infinitely long 
duration, which leads to zero power output. The whole 
cycle is then repeated. In this way heat q from the hot 
bath is converted into work w with stochastic efficiency: 


The macroscopic efficiency is bounded by: 


(w) , Tc 

V= tT - = ^- 7^- 

9 ) Th 


(7) 


This bound is saturated only in the quasi-static regime 
M —>■ 0, with the additional limit x = 1/L ^ 0. For large 
u on the other hand, the set-up will no longer operate 
as a heat engine, since the average amount of delivered 
work becomes negative. The time-reversed process cor¬ 
responds to an expansion of the piston, while being in 
contact with the cold reservoir, followed by a compres¬ 
sion while in contact with the hot reservoir. Under this 
protocol, the system operates as a refrigerator. 


III. QUASI-STATIC LIMIT 

We first focus on the quasi-static limit of the consti¬ 
tutive process, and consider the work delivered by the 
particle, initially at equilibrium and in presence of the 
thermal wall at temperature T, on a piston, moving with 
dimensionless speed u from the initial x = //L to fi¬ 
nal (reduced) location 1. The quasi-static regime corre¬ 
sponds to the limit in which the velocity goes to zero, 
M —?► 0. At lowest order, one obtains the familiar result 
that the work converges to the decrease AF in free en¬ 
ergy. An interesting feature appears when considering 
the first correction, which turns out to be of order u\nu. 


In appendix A, it is shown that the probability distri¬ 
bution Px,u{w) for the work w, upon moving the piston 
from a: to 1 at constant velocity m, is, at this order, given 
by a Gaussian distribution, cf. figure O 

Px.«(w) = -e +0{u). (8) 

's/a'K(7x^u 

The convergence to the Gaussian entails that all cumu- 
lants of higher than second order decrease proportional 
to u in the limit u —> 0: 

+ + (9) 

Remarkably, the reduction in average work due to the 
finite speed of the piston only depends on the reduced 
speed: 

i^)x,u = +■ ( 10 ) 

Similar for the dispersion: 

n 

=—2kBTAF\ —u\nu. (11) 

We stress that the ulnu correction stems from the fact 
that the average time for a particle, initialised in the 
equilibrium state, to reach the thermal wall diverges log¬ 
arithmically as M —?► 0. We note that this result differs 
from the linear dependence found in e.g. low dissipative 
systems [dlj. Furthermore, the convergence to a Gaus¬ 
sian distribution is rather slow, as the non-Gaussian cor¬ 
rections of order u only become negligible compared to 
rtlnw for very small values of u. 

The corresponding result for the time-reversed process, 
i.e. compression of the piston from position 1 to position 
X with dimensionless velocity it, is obtained by replacing 
wlnw with —ulnit and AF by — AF. Hence the corre¬ 
sponding probability distribution, denoted by px,u[w), is 
also Gaussian with first two central moments 


= (^)x.u-2AF 

^X,U — ^X,U‘ 


( 12 ) 

(13) 


We are thus in a position to perform an independent 
check of the above finding, namely by verifying the 
Grooks relation [T^. Since 


Px,u{'>J^) 

Px,u{-w) 


= exp 


= exp 




Mx.« + Mx.u, 

- ^5 - iw - AF) 


we find that the Crooks relation implies: 

'\W)x,u + '\^)x,u ^ _ 


(14) 


(15) 


which is indeed in agreement with Eqs. (nni, (ED and 

ED- 
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A. The Szilard engine 


For the Szilard engine, the distribution of information 
gain and work during one iteration can be written as: 

Pi{Ai,w) = xpx,u{w)SiAi + hix) 

+ (1 - x)pi-x,u{w)S{Ai + ln(l - x)). (16) 

The probability distribution of the stochastic efficiency 
of one iteration can be written as: 


Pi{r]) = xpx,u{-kBTr]\nx) 

+ {1 - x)pi-x,ui-kBTriln{l - x)). (17) 


Combining with Eq. we conclude that the probability 
distribution of the efficiency of the Szilard engine in the 
quasi-static limit can be written as a sum of two Gaussian 
distributions. Furthermore, the macroscopic efficiency is 
found to be, cf. Eq. (ITUl) : 


V 


kBT (In a; + ln(l — a;)) 


1 -I- W —itlnu -I- 0{u). 

TT 


(18) 


Note that the first order correction to the macroscopic 
efficiency is independent of the initial position of the pis¬ 
ton, which leads to universal behaviour in the efficiency 
near the quasi-static limit, as can be seen in figure |4^. 


B. The simplified Carnot heat engine 


The above calculations can also be repeated for the 
heat engine. The joint probability distribution of pro¬ 
duced work and absorbed heat from the hot reservoir in 
one cycle can be written as: 

Piiw,q) = 

/ dwoPx^{wo)p‘' - Wo) f{q - Wo). (19) 

The index h or c appearing in the probability distribu¬ 
tion of the constitutive process refers to the temperature 
of the thermal wall (Th and Tc, respectively) during that 
stage of the cycle. The extra factor y^Th/Tc in the sub¬ 
script of accounts for the fact that the velocity is al¬ 
ways relative to the thermal velocity of the wall, and the 
factor f{q — Wo) accounts for the change in kinetic en¬ 
ergy upon switching the temperature of the heat bath at 
the end of stages 2 and 5. The corresponding probability 
distribution f{q) is given by: 


/(?) = 


27rT, 


S [q- ki 


Thvl - Tcvl 


TvkBVnr, vUsTc ^ fcsTj 2 j ’ 


P 0 . 5 ,n{w/AF) 



Figure 3: (Color online) Gaussian approximation (solid lines) 
of the work distribution for an expansion from I = 0.5 to L = 1 
near the quasi-static limit, together with numerical results, for 
u = 10“® (black crosses), 10“® (orange dots) and 10“^ (purple 
squares). Inset: second (solid line), third (dashed line) and 
fourth (dotted line) cumulant of the scaled work distribution. 


with Ko the modified Bessel function of the second kind 
( 4 ^ . Although there appears to be no simple expression 
for the probability distribution of the stochastic efficiency 
r] = w/q, we can calculate analytically the macroscopic 
efficiency in the quasi-static limit, cf. appendix A: 

_ _ 2{Th - Tc) In a: 

^ 2Th \nx + Tc-Th 

2\/2\nx {2y/T^\nx+ Tc-Th){Th + VT^) 

+ -7^ - (2nin,+T.,-nr -“ “ 

+ 0{u). (21) 

This is confirmed by figure HJd. As mentioned earlier, 
Carnot efficiency can be reached in the limit a: —>• 0: 

^ = 1 - ^ + + V^) “Inu-bO(u), (22) 

see also figure [Ua. 


IV. EFFICIENCY FLUCTUATIONS 

We now turn to the efficiency fluctuations of the en¬ 
gines. The focus here is mainly on the large deviation 
function (LDF) J{r]) of the pro bability distribution Pn{v) 
of the stochastic efficiency : 

J{n) = - lim - lnPn{r]), (23) 

n—^oo Ti 

with n the number of cycles of the engine. For systems 
under time-asymmetric driving, the LDF of the time- 
forward process J{r]) and the LDF of the time-reversed 
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a) \^-v\ b) \riq.s.-V 




Figure 4: (Color online) Macroscopic efficiency of a) the Szi¬ 
lard engine near the quasi-static limit, with x = 0.7 (light 
blue), X = 0.5 (orange) and the time-inverse process of a; = 0.5 
(purple), and b) the heat engine near the quasi-static limit 
with X = 0.5 for Th = l.lTc (red) and Th = 2Tc (green). 


J{v) 



Figure 5: (Color online) — (1/n) lnP„(77) for 10 (blue dotted 
line), 20 (red dashed line) and 50 (dashed dotted green line) 
cycles of the Szilard engine and its time-inverse, with n = 0.1 
and X = 0.7. The purple full curve is the extrapolation to the 
LDF. The macroscopic efficiency is fi = 0.80. The inset shows 

Pioiv)- 

process J ( 77 ) intersect at reversible efficiency ijr [13, : 

Jir^r) = J{-nr)- (24) 

In our numerical simulations, it is of course impossible 
to take the limit of an infinite number of cycles n —> 00 . 
However the convergence to this asymptotic limit is in it¬ 
self an interesting feature. We include the estimate of the 


large deviation function using an extrapolation scheme, 
which is presented in detail in [mil , and which leads 
to accurate results. 

The approach to the LDF for the Szilard engine are 
shown in figure [5l including the result of the aforemen¬ 
tioned extrapolation scheme. The obtained estimate 
of the efficiency LDF reproduces the expected results: 
J{ri) > 0 and reaches zero at macroscopic efficiency; the 
LDF of forward and the reversed process intersect at 
reversible efficiency rjr- Note however that the present 
model differs from the theory discussed in [ 13 , EM Elj , 
because the information gain/loss in the engine cannot 
become zero, even at a stochastic level. Consequently, 
the LDF has no finite plateau at infinity and no power- 
law behaviour is found in the corresponding probability 
distribution. 

The approach to the LDF of the time-forward and 
the time-reversed heat engine are plotted in figure IH] 
The LDFs intersect at Carnot efficiency, as expected. 
The other predictions of the general analysis [l^ are 
reproduced: the maxima of both LDFs have the same 
height J{rj) converges to the same finite plateau value 
for 77 —> ± 00 , and similarly so for the time-reversed pro¬ 
cess. Concerning the finite time regime, we note that the 
crossing of the probability distributions is actually always 
quite close to the reversible efficiency and seems to con¬ 
verge to Carnot efficiency as a power law in function of 
the number of cycles, cf. inset of figure [51 

In order to derive analytical results for the efficiency 
fluctuations, we introduce a slightly modified version of 
the heat engine, and take the limit Tc —>■ 0. When the 
temperature of the thermal wall is zero, the particle sticks 
to it upon collision. The modification is that the hot 
reservoir at temperature Th is present only during an in¬ 
finitesimal short time at the start of the cycle. In this 
way, the particle is launched with a velocity v drawn 
from the Rayleigh distribution (Eq. ([I])) at temperature 
Th- Immediately following the launch, the temperature 
of the thermal wall is dropped to zero. Keeping the relax¬ 
ation step after the expansion of the piston, the particle 
eventually returns to the thermal wall and the piston can 
be returned to its initial position without delivering any 
work. It is clear that the relaxation step at the end of the 
cycle becomes irrelevant, as the particle is immobilised, 
and can be discarded. With this time-symmetric proto¬ 
col, the particle will collide exactly once with the piston 
during each cycle. Furthermore, setting x = 0.5, the 
particle delivers work only when v > 2u. The efficiency 
distribution during one cycle is given by: 
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Pi{r]) = J dv (j){v) ^6{2u — v)S(t]) + 9{v — 2u)S ^ 


( 


= 1-e 


fiv) - V 


-0{1-V^), 


(25) 



J{v) 



Figure 6 : (Color online) — (1/n) lnPn(? 7 ) for 5 (blue dotted 
line), 10 (red dashed line) and 20 (green dashed, dotted line) 
cycles of the heat engine and its time-inverse, with Th = 2 Tc 
(i.e. r}c = 1/2), u = 0.3 and x = 0.5. The purple full curve 
is the extrapolation to the LDF. The macroscopic efficiency 
is given by i; = —0.02 Inset: convergence of the intersections 
efficiency r]* of forward and time-reverse curves to rjc as the 
number n of cycles increases. The dashed line is a power law 
fit of the form ajnP, with a = 5.49 • 10“® and /3 = 0.13 


with f{ri) = 1 -b — rj and 6{x) the Heaviside func¬ 
tion. The LDF can be found by numerically contracting 
the cumulant generating function, as is done in the ap¬ 
pendix. The results are shown in figure [T] It is clear that, 
even at a stochastic level, only efficiencies between 0 and 
1(= rjc) can be reached. Furthermore, the LDF diverges 
at Carnot efficiency and the results from the extrapola¬ 
tion are in good agreement with the analytical formula 
derived in appendix B. 


V. DISCUSSION 


We have analysed the performance of two typical sin¬ 
gle particle engines, the Szilard engine and a simplified 
heat engine. Their analysis is based on the dynamics of 
a single particle moving back and forth between a ther- 
malising wall and a moving piston. For both engines, 
we have shown that the distribution of delivered work 
during one cycle becomes Gaussian in the quasi-static 
limit. More precisely, the corrections of the first two cu- 


Figure 7: (Color online) — (1/n) lnP„( 77 ) for 1 (blue dotted 
line), 2 (red dashed line) and 5 (green dashed, dotted line) 
cycles of the zero-temperature heat engine, with u = 0 . 1 . 
The purple full curve is the extrapolation to the LDF and the 
black curve is the analytical LDF. The macroscopic efficiency 
is given by 77 = 0.23. Inset: Pi{rj). 


mulants in the quasi-static limit scale as it In it, whereas 
the higher order cumulants are proportional to u. For 
the heat engine, both work w and heat q can take on any 
real value, and the fluctuations of efficiency rj reproduce 
the familiar universal features. For the Szilard engine 
however, the information gain per cycle can only take on 
two (positive) values — In a; or — ln(l — x). As a result, 
the asymptotes of the LDF’s for the efficiency diverge, 
since extreme (positive or negative) efficiencies can only 
arise from a corresponding extreme value of the deliv¬ 
ered work. On the other hand, the universal crossing of 
forward and time-reverse LDF’s at reversible efficiency 
1 is reproduced. Our results can easily be extended to 
other systems such as the full Carnot cycle [3l| and the 
Andersen cycle 3^. A question for further research is 
to clarify how the sub-linear u In u behavior fits into the 
linear thermodynamics description of piston engines . 
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Appendix A: Quasi-static limit 


In , it is shown that the probability distri¬ 
bution Px,u{w,t) of the produced work, with r = 

In -I- ut/Lo y/ksT/m ^, Lq being the initial position 
of the piston, t the time, T the temperature, and m the 
mass of the particle, is given by: 


A) = --- -b 0 (m), 

1 - ’tpu{u},X) 


(Al) 


with Px,u{pJ,X) the probability distribution of the deliv¬ 
ered amount of work, with w Fourier transformed to lo 
and T Laplace transformed to A. Furthermore, 

Pi-x,u{^\ A) = 1 - A + iuj -b 0(u), (A2) 

= (A3) 


P f -x ,u{^ 1 ■^) — 

i (l - 540, A) + A) - Mco; A)) , 

(A4) 


( u\^ 

gu{u},X) = / dv(l){v)[l -e 

Jmax(M,0) ^ 

= 1 — uJ —A — iuJ— ujUbT 


—iLoksTuv 


A(A-1)^2 


In u-b O(rt^), (A5) 


with (tq)^ ^ the average modified time before the first 
collision with the thermal wall, {wq)^ ^ the average de¬ 
livered amount of work after the first collision with the 
thermal wall, and (/)(u) the Rayleigh distribution. Filling 
in the last equation in the upper two, gives: 



Figure 8: 3rd cumulant of the work distribution in function 
of u for X = 0.5, with u and ulnu fit. 



Figure 9: 4th cumulant of the work distribution in function 
of u for X = 0.5, with u and ulnw £t. 


so that: 


Pi-,x,u{^', A) = 1 + X\j —uhiu + 0{u). (AlO) 


A) = 1 — V^u (A + iksToj) 

— 2A(A — l)it^ Inw-b O(u^), (A6) 

and 

Pf-x,u{ijJ', X) = + 2(A — 1)m^ Inrt -b O(it^). (A7) 

If the system is initially in its steady state, we have: 

=-Y^Mlnu + 0(u), (A8) 

(A9) 


Combining these results, leads to: 


Px^uip^l A) — 


1 


-b 


A + iiuksT 
2 2iujXkBT -b A^ — iwkBT 

TT (A -b iojkBTY 


u\tlu + 0{u). (All) 


The probability distribution of the produced amount of 
work w and the information consumption Ai can than be 
directly calculated by inverse Laplace and Fourier trans¬ 
forming. Furthermore, Eq. (jA8l) shows that the aver¬ 
age time before the first collision diverges logarithmically 
when u —>■ 0. 


'^'>^o)x,u = 0 ( 4 , 









To numerically check that the higher order cumulants 
of the work distribution go as u instead of ulnu, we 
plot the third and fourth cumulant of the distribution 
of w/AF in function of u and fit to both u and itlnu. 
These plots clearly show that those cumulants are pro¬ 
portional to u, as was discussed above. 

It is also possible to calculate the efficiency of the heat 
engine. The delivered amount of work can be directly 
found from the above calculation: 

iw)x,u = -kB \nx{Th-Tc) 

In a; (Xh + \/ThT^ itlnu -|- 0{u). (A12) 



after which Eq. m is found. 


Appendix B: Zero temperature engine 


The calculation of the cumulant generating function 
or large deviation function of the efficiency of the zero 
temperature engine is a non-trivial problem as efficiency 
is a non-additive quantity. Instead, we shall first turn to 
the cumulant generating function of the work and heat: 


The dissipated heat on the other hand can be written as 
the sum of the delivered amount of work and the heat 
dissipation when the temperature of the heat bath is 
changed: 


A(A, /r) = In = In [A(A, m) + i?(A, m)] , (Bl) 

with: 


{9)x,u = T.) 


— ksThlnx \l + \ —uhiu\ + 0{u), (A13) 


A(A, /r) = ^ (l - exp {2v?{kBTiJ, - l)) , (B2) 


J 




{l-kBT^l)^ 


^-2{l+kBTX)u^ ^ 


“P ““) P* (Vl- W*' -kBT{\ + M)). 


(B3) 


This can easily be numerically contracted to find the LDF 
of the efficiency: 

J{vi) = (B4) 

In particular, we see that in the limit 77 —>■ 1: 

J(77)~ln(l-77). (B5) 
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